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Abstract. Starting from a product initial state, equal-time correlations in 
nonrelativistic quantum lattice models propagate within a lightcone-like causal region. 
The presence of entanglement in the initial state can modify this behaviour, enhancing 
and accelerating the growth of correlations. In this paper we give a quantitative 
description, in the form of Lieb-Robinson-type bounds on equal-time correlation 
functions, of the interplay of dynamics vs. initial entanglement in quantum lattice 
models out of equilibrium. The bounds are tested against model calculations, 
and applications to quantum quenches, quantum channels, and Kondo physics are 
discussed. 


1. Introduction 

Correlations are a quantity of great importance in statistical and condensed matter 
physics, and they represent an essential resource in quantum information science. In 
traditional condensed matter systems, correlations are often measured in scattering 
experiments. More recently, technological advances have established trapped ultracold 
atoms and ions as versatile experimental platforms for the study of many-body quantum 
systems. Owing to the high level of precision and control in such experiments, equal-time 
correlation functions can be measured at atomic spatial resolution and, simultaneously, 
with a temporal resolution much higher than the intrinsic dynamical time scales [1, 2, 3]. 

Theoretically, the creation and propagation of correlations is well understood in the 
case of uncorrelated initial states, and also for exponentially clustered ones. Making use 
of Lieb-Robinson bounds [4], rigorous estimates of the spatial and temporal behaviour 
of equal-time correlation functions have been derived for short-range interacting systems 
with exponentially clustered initial states [5], and also for rather general types of 
interactions and uncorrelated initial states [6]. The picture that emerges from these 
results is that of quasilocality. In the case of short-range interactions this means that 
correlations are approximately (up to exponentially small corrections) conhned to a 
causal region, resembling the lightcone of a relativistic theory; see hgure 1 (left) for an 
illustration. 
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Figure 1. Left: Lieb-Robinson bound for equal-time correlation functions between 
sites i and j, plotted as a function of time t and distance S = d{i,j) between sites. 
Outside the light-coloured cone-like region, correlations are exponentially suppressed. 
Right: Exact analytic results for the absolute value of the connected correlations 
(crfcrj)c in an XX spin chain with nearest-neighbour interactions (37), starting 
from the product initial state (2). Quasilocal behaviour, i.e., spreading with only 
exponentially small effects outside a cone-shaped region, is observed. 


As an example consider a spin chain with nearest-neighbour interactions, as 
sketched in hgure 2 (top left). We are interested in the time evolution of connected 
spin-spin correlations 




( 1 ) 


between lattice sites i and j, where af denotes the x-Pauli matrix at site i. Starting 
from a product state 

l^) = li)*(8)lt)fc (2) 

k^i 

where | t')A: denotes an eigenstate of all connected correlations vanish initially. 
Under the time evolution induced by a Hamiltonian with nearest-neighbour interactions, 
correlations build up and spread in a distance-dependent fashion, as illustrated in hgure 
1 (right) and hgure 2 (bottom left). The larger the distance 5 between i and j, the 
longer it takes for correlations between the sites to build up. 

It should not come as a surprise that the presence of entanglement in the initial 
state can modify this picture. In the same way that performing a measurement 
on one constituent of a Bell pair has an instantaneous ehect on its distant partner, 
so can interactions with entangled entities speed up the propagation of correlations. 
An illustration of this kind of behaviour is given in hgure 2 (right). We consider 
again connected correlation functions {afaj)c of a spin chain with nearest-neighbour 
interactions, this time starting from an initial state 


IV’) — (I t)j-ri| i)j-i + I Di-ril I t)fc 


( 3 ) 
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Figure 2. Exact analytic results for connected correlation functions (erf crpc of the 
XX-dnain with nearest-neighbour interactions (37). Left: Starting from a product 
initial state (2), correlations build up after a time t that scales linearly with the 
distance S between lattice sites i and j. Right: For an entangled initial state where 
initial correlations extend over a distance £ = S — 2. In this case, correlations between 
sites i and j are created on a timescale that is independent of the distance 6 between 
the sites. 


that is mostly of a product form, with the exception of sites i + 1 and j — 1 being 
maximally (Bell) entangled. Due to the almost-product structure, correlations between 
sites i and j are initially vanishing. Under time evolution, the propagation of correlations 
is enhanced by the long-distance entanglement present in the initial state. For the 
specihe, and rather artihcially constructed, initial state (3), correlations build up in a 
distance-independent fashion (hgure 2 bottom right). 

In general, and in particular for physically realistic initial states, the creation of 
correlations will be determined by an interplay of dynamical effects due to interactions 
on the one side, and of initial entanglement on the other side. In this paper we develop 
theoretical tools for a quantitative description of the propagation of correlations in the 
presence of initial entanglement. The main result (22)-(24) is an upper bound on the 
connected correlation function, containing a Lieb-Robinson-type contribution capturing 
the dynamics, and a second term that takes into account initial correlations. Depending 
on the amount and shape of initial correlations, the bound is able to capture the extremes 
of product initial states (as in hgure 2 left) on the one side and long-distance entangled 
initial states (as in hgure 2 right) on the other side, as well as the more involved cases 
in between where the interplay of dynamics and initial entanglement leads to nontrivial 
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propagation patterns. 

The results of this paper apply to a broad class of quantum mechanical lattice 
models (detailed in section 2) and arbitrary initial states, and this generality accounts 
for many potential physical applications. In fact, unless specihcally prepared, a 
product initial state should be considered the exception rather than the rule. Physical 
applications in which long-distance correlated initial states play an important role 
include: 

(a) Quenching away from a quantum critical point. A simple, exactly solvable example 
is a spin-1/2 Ising chain in a transverse magnetic held of strength h [7]. At a critical 
value h = he of the held strength the model undergoes a quantum phase transition 
from a ferromagnetic to a paramagnetic phase. Preparing the system at he in the 
ground state, connected correlation functions between spins at sites i and j decay 
like a power law with the distance between the sites. A quench, i.e., a sudden change 
of the Hamiltonian, then triggers a time evolution, and the propagation of a local 
perturbation will be ahected by the presence of long-distance initial correlations. 

(b) Quantum transport and qubit transfer in spin chains with long-distance entangled 
ground states. Examples are dimerised open chains, some of which are known to 
have ground states with long-distance entanglement between the end points of the 
chain [8] . When such a chain is used as a quantum channel, entanglement-enhanced 
propagation is observed [9], reminiscent of the scenario depicted in hgure 2 (right). 

(c) Building-up of a Kondo screening cloud in the vicinity of an impurity spin. In 
the Kondo model, an impurity spin is coupled to a noninteracting Fermi gas. At 
zero temperature the spatial correlations in the Fermi gas decay like a power law 
with the distance [10]. Starting from an initial state with no correlations between 
impurity and Fermi gas, analytical [11] as well as numerical [12] calculations show 
that correlations build up predominantly in a cone-shaped region in space-time, 
but with a slow (power law) spatial decay outside the cone. 

In section 2 a Lieb-Robinson bound in a rather general setting is reviewed. In sections 
3-5 our main result, a bound on the connected correlation function in the presence of 
long-distance initial entanglement is derived. Specihe types of interactions and initial 
entanglement distributions are discussed in sections 6 and 7. We hnd good qualitative 
agreement when comparing our bounds to model calculations, and also discuss the 
bounds in the context of examples (a)-(c). 

2. Lieb-Robinson bounds 

The strategy is to use Lieb-Robinson bounds to construct an upper bound on equal-time 
correlation functions. As a setting we choose the rather general class of quantum lattice 
models of [6] for which Lieb-Robinson bounds, and also bounds on equal-time correlation 
functions in the case of product initial states, have been derived. The following notation 
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and conditions are similar to those in [6], bnt we have simplified the presentation, and 
in particnlar avoided the C*-algebraic langnage nsed in that reference. 

On a graph A we have a finite-dimensional Hilbert space at each vertex i G A. 
On the tensor prodnct space = <S)i^A Hamiltonian 

veA 

is defined, where the interaction ‘h(X) is a bonnded linear operator acting nontrivially 
only on the part of the Hilbert space that is associated with the snbset X C A. We 
denote the time evolntion of a bonnded linear operator A in the Heisenberg picture as 

A{t) = e'^Me-^^h (5) 


To be able to prove a Lieb-Robinson bound, the interactions ‘h(X) need to decay 
with the spatial separation of the lattice sites in X in a suitable way. This is enforced 
by requiring 


$11 


:= sup 

tie A 


E 

X3i,j 


F{d{t,j)) 


< oo, 


( 6 ) 


where F : [0, cx)) —)• (0, oo) is a positive function characterizing the spatial decay of the 
interactions, and d is the graph distance on A. The requirements on the spatial decay 
function are that F is uniformly summable over A, 


||F|| := sup^F(d(qj)) < oo, 
and that it satisfies 


isA 

jSA 




tieA 


fceA 


F{d{t,j)) 


(7) 

( 8 ) 


On a regular lattice like A = equations (7) and (8) are satisfied for example by 
F{x) oc (1 -h for a > D, which is a suitable choice for pair interactions decaying 
asymptotically for large distances according to a power law with exponent —a. For 
interactions of finite range (like nearest-neighbour interactions), or for exponentially 
decaying interactions, F{x) oc e““'^/(l -f with a > 0 is a suitable choice |. 

Considering bounded linear observables A, B acting nontrivially only on the regions 
X, y C A, respectively, a Lieb-Robinson bound 

|||A(i),S]|| < (9) 

iex jeY 


holds for any f G IR, where 


9{t) 


exp(2||$||C|f|)-l 

exp(2||$||C|f|) 


for d{X,Y) > 0, 
else. 


( 10 ) 


For a proof, see section 2.1 of [6]. 

f As pointed out in [6], a pure exponential F{x) = e~°‘^ does not satisfy (8). 
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The bound (9) can be used to quantify the difference between an observable A(t) that 
is time-evolved with the full Hamiltonian H, and an observable A'{t) obtained by time- 
evolving A with a modihed Hamiltonian H'. To this aim, one can write 


\\A{t)-A'{t)\\ 




'0 


* d 
dr 




-t 


dre-'^" [H - H',A'{T)]e 


/ /iZ/^M AHt 


r\A 

< / dT||l/l'(T).ff-H']|l. (11) 

Jo 

where unitarity of the time evolution operators and the triangle inequality were used (see 
equation (SIO) of [13] or Lemma 3.3 of [6]). The integrand in the last line of equation 
(11) has the form of the commutator on the left-hand side of the Lieb-Robinson bound 
(9), which we can therefore use to further estimate (11). 


4. Decoupled dynamics 

In section 3 we did not specify the way in which the Hamiltonian is modified, i.e., 
how H and H' are related. Since Lieb-Robinson bounds establish quasilocality of the 
time evolution, they can give a particularly useful estimate when comparing the full 
dynamics under a Hamiltonian H with that of a “decoupled” Hamiltonian H' in which 
all interactions between two spatial regions have been eliminated. Quasilocality then 
suggest considering a ball 

Sx := {j e A : d(j,X) < r} (12) 

around the support X of the observable A as one spatial region, and the outside of the 
ball as the other one, and choosing a ball with radius r large enough that the effect of the 
eliminated interaction terms on the time-evolved observable A is small. The decoupled 
Hamiltonian is then given by 

H'= Y. ( 13 ) 

ZeA:ZnSx=0 or ZnS 5,=0 

where denotes the complement of Sx with respect to A. H — H' then contains all 
terms (and only those) that couple Sx to its complement. 

For these choices of H and if', and using the triangle inequality, we can bound the 
integrand in the last line of (11) by 

|||/l'(r),/f-i/']||< ^ |||/l'(r),4(Z)|||. 

ZeA:ZnSx^0 or zn5^^0 


( 14 ) 
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Applying the Lieb-Robinson bound (9) to each of the terms in the sum, and following 
equations (3.10)-(3.15) of [6], one arrives at 

|||.r(T).g-g'lll < 2g(T)||/t||||.i.|| ^+ ^ ^ F{d(i,j]). (15) 

iex 

Inserting this expression into (11), we obtain a bound on ||A(f) — A'(t)||, i.e., a bound 
on the size of the truncation error when comparing the time evolution of A under the 
decoupled Hamiltonian H' to the full time evolution under H. An analogous result is 
obtained for the time evolution of B by using a ball Sy centred around the support Y 
of R. 

5. Spreading of equal-time correlations 

Our main goal is to estimate equal-time connected correlation functions 

(A(t)R(t)), := {A{t)B{t)) - {A{t)){B{t)), (16) 

where (■) = Tr(-p) denotes the quantum mechanical expectation value with respect to 
some initial state p. The strategy is to express the occurring operators in terms of 
differences A{t) — Aft), whose absolute value can be estimated by equations (9) and 
(11). To achieve this, similar to equation (S9) in [13] we write 

(AR)c = {AB) + {A'B) - (A'B) + {A'B') - {A'B') 

'' -V-' '-V-^ 

= 0 =0 

- {{A - A') + {A')) fB - B') + {B')l 

= ((A - A')B) + {A'{B - B')) -{A- A'){B') 

- {A){B - B') + {A'B’) - {A'){B'), (17) 

where the time dependences of the operators have been suppressed. A bound on the 
absolute value of the connected correlator is then given by 

\{AB)f < 2||A - A'||||R|| + 2||A||||R - R'|| + |(A'R'),|. (18) 

In the case of a product initial state, |{A'R')c| is zero as long as the radii r of the 
balls Sx and Sy centred around X and Y, respectively, are non-overlapping, i.e., for 
r < d{X,Y). In this case one recovers the result of [6]. A related result in [5] permits 
exponentially (in space) decaying initial correlations, showing that, in the case of finite- 
range interactions, correlations are restricted to a cone-like region in space-time, with 
only exponentially small corrections outside the causal cone. 

Here we want to allow for arbitrary initial correlations and investigate their effect 
on the creation of correlations in time. The fact that A' and B' evolve under decoupled 
dynamics imposes a restriction on the size of their correlations, which can be seen as 
follows. Divide the total Hilbert space into three factors, 

JYa = M'sx ® ’^Sy ® ’^K\{Sxr\SY)i 


(19) 
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corresponding to different parts of the lattice as indicated by their indices. Since the 
decoupled dynamics does not mix between the factors of this tensor product, we can 
write 

|(A'5')c| = |(A' (g) 5' (g) 1) - (A' (g) 1 (g) 1) (5' (g) 1 (g) 1)1 

<||gl|||| 5 ||Cor(^x:^y), (20) 

with 


Cor(S'x : Sy) := max |(OsxC'sy)c| , 

l|OSxl|,||OSy||<l 


( 21 ) 


where Osx and Osy are observables supported on Sx and Sy, respectively. Cor(S'x : Sy) 
quantifies, for a given state p with respect to which the expectation value on the right- 
hand side of (21) is taken, the amount of correlations between the two regions Sx and 


Combining all the above results we obtain 
\{A{t)B{t)), 


B\ 


< S§rit) := Cor(^x(r) : ^y(r)) 


E E +E E w,j)) 


( 22 ) 


i&Y j&s^{r). 


with 

G(t) (23) 

where the integration of g [as defined in (10)] is elementary. 

At this point we have made explicit the dependence of the right-hand side of (22) 
on the radius r < d{X, Y)/2 of the balls Sx and Sy §. Choosing r small will in general 
reduce the contribution Cor(S'x(r) : S'y(r)) stemming from the initial correlations, but 
will lead to a larger contribution from the Lieb-Robinson term, and vice versa. It is the 
interplay of these two contributions that can lead to interesting propagation patterns 
going beyond those that emerge from product initial states. An optimized bound ^ can 
be obtained by considering r to be t-dependent, and minimizing the right-hand side of 
(22) over r{t) separately for each time t, 


Pillion 


< S^(t) := min 

0<r(t)<d(X,y)/2 


(24) 


This amounts, at any fixed t, to making r{t) just large enough to encompass the causal 
region to which the propagation is essentially restricted, but no larger, in order to reduce 
the contribution from correlations of the initial state. Initial correlations between certain 
regions become relevant only once those regions have been “reached” by the quasilocal 
dynamics. Equation (24), along with (21)-(23), is the main result of this paper. 


§ In principle, different radii could be chosen for Sx and Sy. 
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A better intuition of the implications of the bound (22) can be obtained by specializing 
the result to correlations between single-site observables (e.g., Pauli operators a* in case 
of a spin-1/2 lattice model), and to specific types of interactions (e.g., nearest-neighbour 
or power law decaying interactions). 

Assuming single-site observables Aj and Bj supported at lattice sites i j, (22) 
simplifies to 


Pillion 


< Cor(^,(r) : ^.(r)) 


+ 4G(t) 


F{d{z,k))+ F{d{j,k)) 


fceSf(r) 


fceSj(r) 


(25) 


Estimates of the remaining summations in (25) can be obtained by integral 
approximation. The bounds in the remainder of this section are less tight due to these 
further approximations, but their functional form becomes more evident. 


6.1. Finite-range or exponentially decaying interactions 


As mentioned in section 2, for interactions of finite range (like nearest-neighbour 
interactions) or for exponentially decaying interactions, the function F{x) oc e““^/(l -|- 
x)^+^ with a > 0 is a suitable choice satisfying (7) and (8). By integral approximation 
we can then bound 


E F{d(i,k)) 

k&S^(r) 


- E 

k^A:d{i,k)>r 


^—ad{i.,k) 

(1 d{i, k))^+^ 



(26) 


with a D-dependent constant c > 0. Simplifying also the time-dependence in (25) by 
estimating 


G(t) 


/■'''dr - 1) < 

h Jo 


C+\\F\\ 

2C2 


e2||<i>||c|t|^ 


(27) 


we obtain 


\\A\\\\B\\ 


<CoT{Si{r)-.S,{r)) + 


4c(C'+ ||F||)e2|l^ll^l‘l-' 


(28) 


For the case of an uncorrelated initial state, r = d{i,j)/2 is the optimal choice for 
minimizing the second term on the left-hand side of (28), and one can further estimate 


\\A\\\\B\\ 


< e“(di|-AM)) X const.. 


(29) 


similar to equation (3.1) of [6]. An analogous result, only with different constants v 
and a, holds also for initial states with exponentially clustered correlations [5]. In the 
presence of longer-ranged initial correlations, however, one would need to minimize (28) 
over r in order to obtain a tighter bound. 
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For power law decaying interactions we can choose F{x) oc (1 + x)““ with a > 0. By 
integral approximation one can then bound 

1 . c 


F{d{i,k)) oc 


(1 + d{i, k)) 


a. — j.O'—D 


k£Sf{r) k£A:d{i,k)>r 

with a H-dependent constant c > 0. Inserting (27) and (30) into (25) we obtain 


\{A{t)Bit)), 


<CoT{S,{r):S,{r)) + 


4c(C'+ ||F||)e2|l^ll^l‘l 


(30) 


(31) 


|5|| ---v~.v w~.vw/. 

As in section 6.1, r = d{i,j)/2 is the optimal choice for minimizing the left-hand side 
of (31) in the case of an uncorrelated initial state, yielding 




w|t| 


I rill X const. (32) 

\B\\ - d{t,jy-D ^ ^ 

For power law interactions with a > 2D there exists a sharper Lieb-Robinson-type 
bound due to Foss-Feig et al [14] , which could be used to derive a bound on equal-time 
correlation functions along the same lines as above. 


7. Examples of bounds for long-distance correlated initial states 

One-dimensional examples are presented, illustrating how correlations of the initial state 
modify the creation and propagation of correlations between initially uncorrelated lattice 
sites. 


7.1. Connecting with an entangled pair 
Here we consider an initial state 

iV’) = (it)-fcU)fc + i;)-fcit)fc) (g) It)™ (33) 

m^—k,k 

that is mostly of product form, except for a maximally entangled pair at sites —k and k. 
We are interested in the time evolution of the connected correlation function {afpjDc 
between lattice sites —i and i that are a distance 6 = 2i apart. For the initial state (33) 
we have 

Cor(S'_i(r) : Si{r)) = {(po\a^_j^al\(po)^ 0(r - \i - k\) = 0(r - \i - k\) (34) 

for r < i, where 0 denotes the Heaviside step function. Assuming a lattice model with 
nearest-neighbour interactions, we combine (28), (29), and (34) to obtain 

\{(rP{t)a^{t))c\ < min|l,min(0(r-|i-A;|)-l-ce“(’^l'l-’'))} 

= min{l,ce“(^l*l-l*-^l)} , (35) 

where we have also included the bound 
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Figure 3. Left: Contour plot of the bound (35) for the initial state (33) and parameter 
values 2k = 10 and a = v = c = 1. Right: Exact analytic results for the absolute value 
of equal-time correlations in an XX-chain with nearest-neighbour interactions (37), 
starting from the initial state (33). In both plots, cone-like spreading of correlations 
similar to the behaviour in the absence of initial correlations is observed, but with an 
offset (i.e., a shift to higher values of S) by 2k. 


which follows from the Cauchy-Schwarz inequality. The contour plots in hgure 3 (left) 
illustrate the creation of correlations in time and as a function of the distance 6 between 
spins. The correlated pair leads to an effective reduction of the distance d{i,j) between 
the spins by the distance 2k between the correlated sites: Mediated by the entanglement 
of the initial state, spins that are 2k 1 sites apart “feel” each other as if they were 
neighbours. Accordingly, the time to transmit a signal across such an entangled quantum 
channel is reduced by 2/c/n, where v is the Lieb-Robinson velocity occurring in (35). 

To assess how well the bound compares to the actual dynamics, we investigate 
the time-evolution of the state (33) under an XX-Hamiltonian with nearest-neighbour 
interactions, 

( 37 ) 

i 

where we set J = 1. As shown in figure 3 (right), correlations spread in the interior 
of a cone, with a spatial offset of 2k = 10 compared to the case of an uncorrelated 
initial state. This confirms that the bound (35) reproduces the qualitative features 
of the dynamics, although, as expected, the velocity at which magnon quasiparticles 
propagate in the XX-chain is slower than the estimated velocity v occurring in the 
estimate (35). The bound and the model calculation shown in hgure 3 can be seen as 
simplifted illustrations of the enhancement of quantum transport and qubit transfer in a 
long-distance entangled quantum channel, as mentioned in item (b) of the Introduction. 
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Figure 4. Top: Sketch of a spin chain with power law clustering of correlations (38). 
Centre left: Contour plot of the logarithm of the bound (39) for power law-correlated 
initial states, with parameters x = 2, o = u = ci = C 2 = 1. Bottom left: As above, 
but for fixed values of t and in a log-log representation. Enter right: Exact analytic 
results for the logarithm of the absolute value of the equal-time correlation function 
for a Kondo impurity coupled to a Fermi sea at zero temperature. Bottom right: As 
above, but for fixed values of t and in a log-log representation. All these plots show 
that correlations spread inside a cone, where the spatial decay outside the cone follows 
a power law, as is visible in the log-log plots. 


1.2. Power law clustering of initial correlations 

Items (a) and (c) of the Introduction describe two possible scenarios of physical interest 
where initial states with power law-clustered correlations arise (illustrated in hgure 4 
top). In that case the connected correlations initially satisfy 

|{<y>c| < 




(38) 
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with some exponent x ^ 0- Assuming again a chain of spin-1/2 degrees of freedom with 
nearest-neighbour interactions, we combine (28), (29), and (38) to obtain the bound 

< min 11,mm | . (39) 

As illustrated in figure 4 (left), this bound shows cone-like propagation, but, in contrast 
to the case without initial correlations, the spatial decay outside the cone follows a 
power law with exponent y (instead of an exponential decay). 

As a physical illustration of this kind of propagation behaviour, we borrow results 
from Medvedyeva et al [11] on the spatiotemporal build-up of the Kondo screening cloud. 
The authors of that paper study correlations between the spin of an impurity and the 
spin of a conduction electron in the three-dimensional Kondo model. Strictly speaking 
this model does not satisfy the conditions under which the Lieb-Robinson bound (9) has 
been proved, and hence our bound (24) does not apply. Proving Lieb-Robinson bounds 
for general bosonic or fermionic hopping models turns out to be elusive, but model 
calculations indicate that the majority of such models with local interactions nonetheless 
do show lightcone dynamics. Counterexamples exist, but they require careful design [15], 
and we have good reason to believe that the Kondo model can well serve as an example 
illustrating the physics described by the bound (24). 

The initial state | t) ® |FS) used in [11] is a product of the impurity spin 
and the Fermi sea of the conduction electrons. While the impurity spin is initially 
uncorrelated with the conduction electrons, the conduction electrons themselves are 
spatially correlated among one another. At zero temperature these initial correlations 
decay like a power law in space. Analytic expressions are then obtained for the 
correlation functions of the Kondo model at the Toulouse point (i.e., for a special value 
of one of the coupling constants in the Kondo Hamiltonian), see equations (9) and (19)- 
(23) of [11]. Numerically evaluating these equations in the zero-temperature limit, we 
obtain the spatiotemporal spreading of the equal-time correlation functions plotted in 
figure 4 (right). The correlations are sharply peaked on the boundary of the lightcone. 
The spatial decay outside the cone follows a power law, in agreement with the bound 
shown in figure 4 (left). 

8. Conclusions 

In this paper we have developed theoretical tools, applicable to a broad class of 
quantum lattice models, for the description of the propagation of correlations in the 
presence of initial entanglement. The main result (22)-(24) is an upper bound on the 
connected correlation function, containing a Lieb-Robinson-type contribution capturing 
the dynamics, and a second term that takes into account initial correlations. Depending 
on the amount and shape of initial correlations, the bound is able to capture the extremes 
of product initial states (as in figure 2 left) on the one side and long-distance entangled 
initial states (as in figure 2 right) on the other side, as well as the more involved cases 
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in between where the interplay of dynamics and initial entanglement leads to nontrivial 
propagation patterns. 

The essential prerequisite for the proof is that some kind of Lieb-Robinson bound 
can be derived for the system under consideration. We have here used the setting of 
[6] (detailed in section 2), which includes short- as well as long-range interacting lattice 
models, but other settings may be used to either allow for different or more general 
types of interactions, or for sharper bounds [16, 17, 18, 19]. The main idea of the proof 
is to divide the lattice into two regions, one consisting of two disjoint spheres Sx and Sy 
of radius r, centred around the lattice sites i and j for which the connected correlation 
function (ufcrpc is to be estimated; and the other region being the complement of the 
two spheres. The size of the correlations can then be bounded by two contributions: The 
first two terms on the right-hand side of (18) account for the propagation of correlations 
due to interactions, and they decrease with increasing r; the third term on the right- 
hand side of (18) accounts for initial correlations between Sx and S'y, and it increases 
with r. An optimal bound is then found by choosing, for each time t and given sites 
i and j, the optimal value of r for which the bound becomes minimal, as in the final 
result (24). 

Owing to its generality, many potential applications of this result can be envisaged. 
Physical applications in which entangled initial states have a strong effect on the 
propagation of correlations include quenches away from a quantum critical point, 
quantum transport and qubit transfer in spin chains with long-distance entangled ground 
states, or the building-up of a Kondo screening cloud at zero temperature. Some of these 
examples, or simplified toy models of them, have been discussed in section 7. While 
absolute magnitudes and propagation velocities are overestimated (as is generally the 
case when using Lieb-Robinson bounds), the bounds on correlation functions derived in 
this paper show good qualitative agreement with model calculations for entangled initial 
states. These comparisons were made for integrable models where exact analytic results 
are available, but we have no reason to believe that the propagation of correlations 
will in general be substantially different in nonintegrable systems. In certain instances, 
however, localisation effects can strongly suppress the propagation of correlations. In 
these cases the bound (24) should be based on a sub-ballistic Lieb-Robinson bound like 
the one derived for disordered spin systems in [20]. 

We believe that the results here provide useful descriptions of the propagation pat¬ 
terns to be expected in a variety of physical situations of interest, and they should be 
particularly expedient for larger system sizes where numerical simulations are out of 
reach. 
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